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, | ■ Abstract 

\ Theoretical understanding of experimental results from relativistic heavy- 

ion collisions requires a microscopic approach to the behavior of QCD n- 
point functions at finite temperatures, as given by the hierarchy of Dyson- 
Schwinger equations, properly generalized within the Matsubara formal- 
ism. The convergence of sums over Matsubara modes is studied. The 
technical complexity of finite-temperature calculations mandates model- 
ing. We present a model where the QCD interaction in the infrared, 
nonperturbative domain, is modeled by a separable form. Results for the 
mass spectrum of light quark flavors (u, d, s) and for the pseudoscalar 
bound-state amplitudes at finite temperature are presented. 
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C2r ■ 1 Introduction 

O 

r"| ■ While the experiments at RHIC [TJ [5] advanced the empirical knowledge of the 

Oh! hot QCD matter dramatically, the understanding of the state of matter that has 

been formed is still lacking. For example, the STAR collaboration's assessment 
\ [5] of the evidence from RHIC experiments depicts a very intricate, difficult-to- 

' understand picture of the hot QCD matter. Among the issues pointed out as 

important was the need to clarify the role of quark-antiquark (qq) bound states 
continuing existence above the critical temperature T c , as well as the role of the 
chiral phase transition. 

■ Both of these issues are consistently treated within the Dyson-Schwinger 

(DS) approach to quark-hadron physics. Dynamical chiral symmetry breaking 
(DChSB) as the crucial low-energy QCD phenomenon is well-understood in 
the rainbow-ladder approximation (RLA), a symmetry preserving truncation 
of the hierarchy of DS equations. Thanks to this, the behavior of the pion 
mass is in accord with the Goldstone theorem: the pion mass shows correct 
behavior while approaching the chiral limit, as seen on Fig. [T] This correct 
chiral behavior is a general feature of the DS approach in RLA, and not a 
consequence of our specific model choice. For recent reviews of the DS approach, 
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Figure 1: Correct chiral behavior of the pion mass close to chiral limit: M% oc 
fh q , where fh q is the bare quark mass parameter. It is akin to the notion of the 
current quark mass in QCD and gives the extent of the explicit chiral symmetry 
breaking as opposed to DChSB - see Eq. ([I]) below. 



see, e.g., Refs. [3j|4], of which the first [3] also reviews the studies of QCD DS 
equations at finite temperature, started in [5]. Unfortunately, the extension of 
DS calculations to non- vanishing temperatures is technically quite difficult. This 
usage of separable model interactions greatly simplifies DS calculations at finite 
temperatures, while yielding equivalent results on a given level of truncation 
O [7] . A recent update of this covariant separable approach with application to 
the scalar a meson at finite temperature can be found in [8]. Here, we present 
results for the quark mass spectrum at zero and finite temperature, extending 
previous work by including the strange flavor. 



2 The separable model at zero temperature 

The dressed quark propagator S q (p) is the solution of its DS equation [3111], 

S.ip)- 1 = < 7 -p + m q + | J -0^ g 2 Df y (p - t)^S q {l) lv , (1) 

while the qq 1 meson Bethe-Salpeter (BS) bound-state vertex T qq ' (p, P) is the 
solution of the BS equation (BSE) 

-X(P 2 )T qql (p,P) = ^g 2 Df u (p~e) ltl S q (£ + )r qql (£ 1 P)S q (e^ 1 (2) 

where Df u (p-l) is an effective gluon propagator modeling the nonperturbative 
QCD effects, rh q is the current quark mass, the index q (or q') stands for the 
quark flavor (u,d or s), P is the total momentum, and l± — I ± P/2. The 
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chiral limit is obtained by setting rh q = 0. The meson mass is identified from 
X(P 2 = —A/ 2 ) = 1. Equations (JTJ) and ([2]) are written in the Euclidean space, 
and in the consistent rainbow-ladder truncation. 
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Figure 2: The sum in Eq. as the function of number of the Matsubara modes 
included in summation at temperature T = 1 MeV. Eq. |[8]) is normalized to 
value calculated with enough Matsubara modes (n = 1000) to achieve prescribed 
numerical precision. 

The simplest separable Ansatz which reproduces in RLA a nonperturba- 
tive solution of ([1]) for any effective gluon propagator in a Feynman-like gauge 
g 2 Df v (p - £) -> d^D(p 2 , £ 2 ,p ■ £) is C| 



D(p 2 ,t 2 ,p-l) = D F (p 2 )F (t 2 ) + D 1 F 1 (p 2 )(p-£)F 1 (£ 2 ) . 



(3) 



This is a rank-2 separable interaction with two strength parameters Di and 



corresponding form factors Ti(p ), ■ 

constrained to the solution of the DSE for the quark propagator ^ 
Sgip)- 1 = il ■ P A q (p 2 ) + B q (p 2 ) = Z-\p 2 )[i~j -p + mgip 2 )}, 



1,2. The choice for these quantities is 

(4) 



where m q (p 2 ) = B q {p 2 ) / A q {p 2 ) is the dynamical mass function and Z q (p 2 ) = 
A~ 1 (p 2 ) the wave function renormalization. Using the separable Ansatz Q in 
(JTJ, the gap equations for the quark amplitudes A q (p 2 ) and B q (p 2 ) read 



B q (p 2 ) - m q 
[A q ( P 2 ) 1] 
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b q T (p 2 ), (5) 
■aMp 2 ). (6) 
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Once the coefficients a q and b q are obtained by solving the gap equations ||SJ) 
and ([6]), the only model parameters remaining are fh q and the parameters of 
the gluon propagator, to be fixed by meson phenomenology. 
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Figure 3: The sum in Eq. ([5]) as the function of number of the Matsubara modes 
included in summation at temperature T — 10 MeV. Eq. ([5]) is normalized to 
value calculated with enough Matsubara modes (n = 100) to achieve prescribed 
numerical precision. 



3 Extension to finite temperature 

The extension of the separable model studies to the finite temperature case, 
T ^ 0, is systematically accomplished by a transcription of the Euclidean quark 
4-momcntum viap — > p n = (u> n ,p), where uj n = (2n-\-\)irT are the discrete Mat- 
subara frequencies. In the Matsubara formalism, the number of coupled equa- 
tions represented by dXJ) and @ scales up with the number of fermion Matsubara 
modes included. For studies near and above the transition, T > 100 MeV, using 
only 10 such modes appears adequate. Nevertheless, the appropriate number 
can be more than 10 3 if the continuity with T = results is to be verified. 
Convergence of the sum in the Eq. ([5]) is shown on Figs. HHD The effective 
qq interaction, defined in the present paper by the Ansatz and the form 
factors (|12p below, will automatically decrease with increasing T without the 
introduction of an explicit T-dependence which would require new parameters. 

The solution of the DS equation for the dressed quark propagator now takes 
the form 

S~ 1 (p» ,T) = ij-pA q (pi , T) + ^ 74Wrl C q (p 2 n ,T) + B q (pi , T) , (7) 
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Figure 4: The sum in the Eq. ([8|) as function of number of the Matsubara modes 
included in summation at temperature T — 100 MeV. Eq. ([8|) is normalized to 
value calculated with enough Matsubara modes (n = 10) to achieve prescribed 
numerical precision. 



where p 2 n — u 2 +p 2 and the quark amplitudes B q (p 2 l ,T) — m q + b q (T)J r o(p 2 l ) , 
A q {pi,T) = 1 + a q (T)ri(jfi), and C q {p 2 n ,T) = 1 + c,(T).Fi(p£) are defined by 
the temperature-dependent coefficients a q (T),b q (T), and c q (T) to be determined 
from the set of three coupled non-linear equations 

a q (T) = f ^ ^i(^)p 2 [l + o,(T)^i(p£)] d- 1 ^, T) , (8) 

n 

C(? (T) = ^T£ / ^^ l( p2 )tJ 2 [1 + C9(Tmb 2 )]d -l b 2 jr)j (9) 
n ^ ' 

W = ^pT£ / ^^)[a,H(W fl K 1 (p!l,r)- (io) 

The function d q (p 2 l ,T) is the denominator of the quark propagator S q (p n ,T), 
and is given by 

d,(p* , r) = P 2 4(Pn> + "lC 2 q (pi, T) + B 2 (p 2 nl T). (11) 

The procedure for solving gap equations for a given temperature T is the 
same as in T = case, but one has to control the appropriate number of 
Matsubara modes as mentioned above. 
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Figure 5: The p 2 dependence of the form factors Tv(p 2 ) and T\{p 2 ). 

4 Confinement and Dynamical Chiral Symme- 
try Breaking 

If there are no poles in the quark propagator S q (p) for real timelike p 2 then 
there is no physical quark mass shell. This entails that quarks cannot prop- 
agate freely, and the description of hadronic processes will not be hindered 
by unphysical quark production thresholds. This sufficient condition is a vi- 
able possibility for realizing quark confinement [7|. A nontrivial solution for 
B q (p 2 ) in the chiral limit (mo = 0) signals DChSB. There is a connection be- 
tween quark confinement realized as the lack of a quark mass shell and the 
existence of a strong quark mass function in the infrared through DChSB. The 
propagator is confining if m 2 (p 2 ) ^ —p 2 for real p 2 , where the quark mass func- 
tion is m q {p 2 ) = B q (p 2 )/A q (p 2 ). In the present separable model, the strength 
b q = B q (0), which is generated by solving Eqs. §5§ and ©, controls both con- 
finement and DChSB. At finite temperature, the strength b q (T) for the quark 
mass function will decrease with T, until the denominator of the quark 
propagator can vanish for some timelike momentum, and the quark can come 
on the free mass shell. The connection between deconfinement and disappear- 
ance of DChSB is thus clear in the DS approach. Also the present model is 
therefore expected to have a deconfinement transition at or a little before the 
chiral restoration transition associated with ba(T) — ► 0. 

The following simple choice of the separable interaction form factors (graph- 
ically represented on Fig. [5]), 

Fo(p ) = cM-P Mo) , MP ) = l + eM{p 2_ p 2 )/Al) , (12) 



G 



is used to obtain numerical solutions which reproduce very well the phenomenol- 
ogy of the light pseudoscalar mesons and generate an acceptable value for the 
chiral condensate. 

The resulting quark propagator is found to be confining and the infrared 
strength and shape of quark amplitudes A q (p 2 ) and B q (p 2 ) are in quantitative 
agreement with the typical DS studies. Gaussian-type form factors are used as 
a minimal way to preserve these properties while realizing that the ultraviolet 
suppression is much stronger than the asymptotic fall off (with logarithmic 
corrections) known from perturbative QCD and numerical studies on the lattice 




Figure 6: The p 2 dependence (at T = 0) of the dynamically generated quark 
masses m s (p 2 ), m u {p 2 ) for, respectively, the strange and the (isosymmetric) non- 
strange case. 



5 Bound-state amplitudes 

With the separable interaction, the allowed form of the solution of Eq. ([2]) for 
the pseudoscalar BS amplitude is [6] 

r P s(£; P) = 75 (iE PS (P 2 ) + fF PS {P 2 )) T {i 2 ). (13) 

The dependence on the relative momentum t is described only by the first form 
factor Tq{1 2 \ The second term T\ of the interaction can contribute to BS 
amplitude only indirectly via the quark propagators. The pseudoscalar BSE 
([2]) becomes a 2 x 2 matrix eigenvalue problem K.(P 2 )V = A(P 2 )V where the 
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eigenvector is V = (Eps,Fps). The kernel is 

K ij (P 2 ) = -^tv s J -^^{^[USS+^SS-)] , (14) 

where i, j = 1,2 denote the components of t — (175, 75 f) and i = (£75, —75 f/2P 2 ). 
The separable model produces the same momentum dependence for both am- 
plitudes (containing Fps and Eps) in the BS amplitude (fTB")) : the dependence 
of the quark amplitude B q (£ 2 ). Goldstone's theorem is preserved by the present 
separable model; in the chiral limit, whenever a nontrivial gap-equation solution 
for B q (p 2 ) exists, there will be a massless pion solution to (TlH) . 

The normalization condition for the pseudoscalar BS amplitude can be ex- 
pressed as 

x S q (e+)^Ps(i;K)S q (£-)]\ p2=K2= _ M2ps . (15) 

Here T(q; K) is the charge conjugate amplitude [C^ 1 T(— q, K)C] 1 , where C = 
7274 and the index t denotes a matrix transpose. Both the number of colors N c 
and light flavors Nf are 3. 

At T = the mass-shell condition for a meson as a qq bound state of the 
BSE is equivalent to the appearance of a pole in the qq scattering amplitude as 
a function of P 2 . At T ^ in the Matsubara formalism, the 0(4) symmetry is 
broken by the heat bath and we have P — > (Q m , P) where fl m = 2mirT. Bound 
states and the poles they generate in propagators may be investigated through 
polarization tensors, correlators or Bethe-Salpeter eigenvalues. This pole struc- 
ture is characterized by information at discrete points Q m on the imaginary 
energy axis and at a continuum of 3-momenta. One may search for poles as a 
function of P 2 thus identifying the so - called spatial or screening masses for 
each Matsubara mode. These serve as one particular characterization of the 
propagator and the T > bound states. 

In the present context, the eigenvalues of the BSE become A(P 2 ) — > A(fi^, P 2 ; T). 
The temporal meson masses identified by zeros of 1 — A(£1 2 ,0;T) will differ in 
general from the spatial masses identified by zeros of 1 — A(0, P 2 ;T). They are 
however identical at T = and an approximate degeneracy can be expected 
to extend over the finite T domain, where the 0(4) symmetry is not strongly 
broken. 

The general form of the finite-T pseudoscalar BS amplitude allowed by the 
separable model is 

r PS (q n ;P m ) = l5 {iE PS {Pl)+^n m Fp S {Pl) 

+ i-PFps{P 2 m )) Mil)- (16) 

The separable BSE becomes a 3 x 3 matrix eigenvalue problem with a kernel 
that is a generalization of Eq. (| 14|) . In the limit Vl m — > 0, as is required for the 
spatial mode of interest here, the amplitude Fps = £l m Fps is trivially zero. 
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6 Results 



Parameters of the model are completely fixed by meson phenomenology calcu- 
lated from the model as discussed in [3 [8]. In the nonstrange sector, we work 
in the isosymmetric limit and adopt bare quark masses fh u = rhd = 5.5 MeV 
and in strange sector we adopt m s = 115 MeV. Then the parameter values 

A = 758 MeV, Ai = 961 MeV, p = 600 MeV, (17) 

L> Ao = 219, DiAf = 40, (18) 

lead, through the gap equation, to a Ut d = 0.672, b u ,d = 660 MeV, a s = 0.657 
and b s = 998 MeV i.e., to the dynamically generated momentum-dependent 
mass functions m q (p 2 ) shown in Fig. In the limit of high p 2 , they converge 
to rh u and fh s . However, at low p 2 , the values of m u (p 2 ) are close to the typical 
constituent quark mass scale ~ m p /2 ~ mjv/3 with the maximum value at p 2 = 
0, m u (0) = 398 MeV. The corresponding value for the strange quark is m s (0) = 
672 MeV. Fig. [5] hence shows that in the domain of low and intermediate 

P 2 ^5 1 GeV 2 , the dynamically generated quark masses are of the order of typical 
constituent quark mass values. 

Thus, the DS approach provides a derivation of the constituent quark model 
[TU] from a more fundamental level, with the correct chiral behavior of QCD. 
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Figure 7: The temperature dependence of B S (Q), B u {0) and Bq(0), the scalar 
propagator functions at p 2 — 0, for the strange, the nonstrange and the chiral- 
limit cases, respectively. The temperature dependence of the chiral quark- 

1 /3 

antiquark condensate, — (qq) , is also shown (by the lowest curve). Both chiral- 

1 /3 

limit quantities, i?o(0) and —(qq) , vanish at the chiral-symmetry restoration 
temperature T ch = 128 MeV. 
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Obtaining such dynamically generated constituent quark masses, as previous 
experience with the DS approach shows (see, e.g., Refs. such as [3j|H[T0]), * s 
essential for reproducing the static and other low-energy properties of hadrons, 
including decays. (We would have to turn to less simplified DS models for 
incorporating the correct perturbative behaviors, including that of the quark 
masses. Such models are amply reviewed or used in, e.g., Refs. [3J 21 [TO1 ITT] , 
but addressing them is beyond the present scope, where the perturbative regime 
is not important.) 

Another important result related to the dynamically generated, dressed 
quark propagator, is the chiral quark-antiquark condensate (qq)o- For the pa- 
rameter values quoted above, we obtain the zero-temperature value (qq)o = 
(— 217 MeV) 3 , which practically coincides with the standard QCD value. 

The extension of these results to finite temperatures is given in Figs. [3 
[51 Very important is the temperature dependence of the chiral-limit quantities 
Bq(0)t and (qq)o(T), whose vanishing with T determines the chiral restoration 
temperature Teh- We find Tch = 128 MeV in the present model. 



m q (0)[GeV] 




Figure 8: The temperature dependence of m s (0), m u (0) and mo(0), the dynam- 
ically generated quark masses at p 2 = for the strange, the nonstrange and the 
chiral-limit cases, respectively. 

The temperature dependences of the functions giving the vector part of 
the quark propagator, A u ^ s (0)t and C u , s (0)t, are depicted in Fig. [9j Their 
difference is a measure of the 0(4) symmetry breaking with the temperature T. 

The results for pseudoscalar Eps and pseudovector Fps amplitudes can be 
seen on Fig.[lOl The pseudovector amplitude for pion is significantly different 
from zero but decreases rapidly above the transition. 

The presented model, when applied in the framework of the Bethe-Salpeter 
approach to mesons as quark-antiquark bound states, produces a very satisfac- 
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Figure 9: The violation of 0(4) symmetry with T is exhibited on the example 
of A U>S {Q,T) and C r u ,,(0,T). 
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Figure 10: The temperature dependence of pseudoscalar covariants in BS am- 
plitude 



tory description of the whole light pseudoscalar nonet, both at zero and finite 
temperatures [12) . The masses and decay constants of the pseudoscalar nonet 
mesons at T = are summarized and compared with experiment in Table [1] 
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The first three rows in Table[JJgive the masses Mps and decay constants fps 
of the pscudoscalar qq' bound states PS = tt + , K + , and ss resulting, through 
the BS equation ||2J|, from the separable interaction ©. The ss pseudoscalar 
meson is a useful theoretical construct, but is not realized physically, at least 
not at T = 0. (It is therefore not associated with any experimental value in 
this table. Also note that the unphysical mass M sS given in this table does not 
include the contribution from the gluon anomaly.) The parameter values [ (|17|) 
and (fT8"|) in the effective interaction D c ^ vl and the bare quark masses rh Ul d = 5.5 
MeV and fh s = 115 MeV] are fixed by fitting the pion and kaon masses and 
decay constants. These masses and decay constants are the input for the de- 
scription of the 77—77' complex [TT] . More precisely, 77 and 77' masses are obtained 
by combining the contributions from the non-Abelian (gluon) axial anomaly 
with the non-anomalous contributions obtained from the results on the masses 
of n,K, and the unphysical ss pseudoscalar [TT]. For this procedure, it is es- 
sential that we have the good chiral behavior of our qq' bound states, which are 
simultaneously also the (almost-) Goldstone bosons, so that 

M^, = Const (m q + rh q >) , (19) 

as seen in Fig. [TJ For example, Eq. (|T9|) guarantees the relation M% + M% = 
2Mj C which is utilized in Refs. [TTJ Q2] in the treatment of the 77-77' complex. 
Indeed, the concrete model results for M v ,Mk and M sS in Table [JJ obey this 
relation up to \%. 



PS 


Mps 




fps 


rCxp 
J PS 


7T+ 


0.140 


0.1396 


0.092 


0.0924 ±0.0003 


K+ 


0.495 


0.4937 


0.110 


0.1130 ±0.0010 


ss 


0.685 




0.119 




77 


0.543 


0.5473 






77' 


0.933 


0.9578 







Table 1: Results on the pseudoscalar mesons at zero temperature, T = 0, and 
comparison with experiment (where appropriate). All results are in GeV. 

Especially interesting is the temperature behavior of the 77-77' complex, where 
the results for the meson masses differ very much for various possible relation- 
ships between the chiral restoration temperature T^h and the temperature of 
melting of the topological susceptibility, denoted by T x . The once favored sce- 
nario of Pisarski and Wilczek [T3], where 77 would smoothly evolve with T to 
purely non-strange ?7jvs, 

\ Vns ) = ^(l««> + \dd)) = ^|r? 8 ) + ^1 1%) , (20) 
and 77' to purely strange 775, 

\Vs) = \ss)=- ] H\ m ) + ^=\vo) , (21) 
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would occur only when T x is significantly below To- Such a case is depicted 
in Fig. [ITJ where T x = 2/3 Tch- In this case, around the chiral restoration 
temperature rj becomes quite light, and one would expect an increase of the 
relative multiplicity of rj mesons around Tch- Nevertheless, the possibility that 
T x < Tch is nowadays disfavored by the lattice results on the temperature 
dependence of the topological susceptibility [13] . 

On the other hand, for T x ~ Tch or T x > Tch we find that 77 never becomes 
light, while rf even becomes very heavy [12) . Thus, for the relationships favored 
by the lattice, our results [T2] indicate so strong suppression of rf around the 
chiral restoration temperature, that it may constitute a useful signal from the 
hot QCD matter. 

M P [GeV], 1>2/3T ch 




0.5 1 1.5 2 



Figure 11: The relative temperature dependence, on T/T x , of the pseudoscalar 
meson masses for T x = 2/3 Tch- The three variously dashed curves represent 
the pseudoscalar meson masses which do not receive contributions from the 
gluon anomaly: M W (T), M K (T) and M sS (T), depicted by the long-dashed, short 
dashed and dash-dotted curves, respectively. The lower solid curve is M V (T), 
and the upper solid curve is M T? /(T). The lower and upper dotted curves are 
the masses of r\ ns and ijs- The thin diagonal line is twice the zeroth Matsubara 
frequency, 2ttT. This is the limit to which meson masses should ultimately 
approach from below at still higher temperatures, where qq states should totally 
dissolve into a gas of weakly interacting quarks and antiquarks. 
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